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ABSTRACT 

-  The  importance  of  appropriate  stochastic  models  in  choosing  efficient 
methods  of  statistical  analysis  is  discussed.  The  fitting  to  data  of 
Seasonal  Autoregressive  Moving  Average  models  is  described  and  it  is  shown 
how  trend  may  be  estimated  in  an  appropriate  class  of  models  of  this  kind. 
The  procedure  is  illustrated  for  a  model  fitted  to  a  money  supply  series 
published  by  the  Federal  Reserve  Board.  Error  limits  are  calculated. 

In  a  series  of  Appendices  the  properties  of  the  adaptive  coefficients 
which  determine  the  trend  estimates  are  derived. 
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SIGNIFICANCE  AND  EXPLANATION 

Seasonal  fluctuations  in  time  series  often  obscure  information  about 
trend.  Thus  it  may  be  known  that  recruitment  for  the  Army  regularly 
decreases  in  certain  months  of  the  year  and  therefore  that  reductions  in 
these  months  are  not  to  be  regarded  as  indicating  any  real  change  in  the 
recruitment  situation.  Historically  problems  of  this  kind  have  been 
dealt  with  by  making  " seasonal  adjustments ”  but  the  manner  in  which  these 
adjustments  were  made  has  been  somewhat  arbitrary  is  occasionally 
misleading. 

It  is  argued  in  this  paper  that  the  problem  can  best  be  dealt  with 
not  by  seasonal  adjustment  at  all  but  by  direct  estimation  of  the  current 
trend  using  an  appropriate  stochastic  model  which  has  been  fitted  to  the 
series.  The  method  has  the  advantage  that,  provided  appropriate  pre¬ 


cautions,  which  are  discussed,  are  taken  to  ensure  that  the  model  is 
adequate,  the  estimates  of  trend  will  have  efficient  properties  and  they 
will  employ  the  data  in  a  manner  most  appropriate  to  the  particular 
series  under  study. 


The  responsibility  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MRC,  and  not  with  the  authors  of  this  report. 


ESTIMATING  CURRENT  TREND  AND  GROWTH 
RATES  IN  SEASONAL  TIME  SERIES 


George  E.P.  Box  and  David  A.  Pierce 

1.  INTRODUCTION 

For  decades  the  seasonal  adjustment  of  time  series  has  been  widely 
practiced  and  even  more  widely  discussed.  The  goal  of  such  adjustment  is 
presumably  to  facilitate  the  elucidation  and  interpretation  of  other  system¬ 
atic  aspects  of  the  series,  i.e.,  of  the  trend.  Thus,  in  one  sense  the 
title  of  this  paper  is  another  name  for  seasonal  adjustment.  However,  our 
approach  differs  in  several  important  respects  from  the  Census  Xll  procedure 
(Shiskin,  Young,  and  Musgrave,  1967)  and  other  traditional  approaches  to 
deseasonalization.  Most  importantly,  it  is  linked  to  building  a  model  for 
the  series,  thereby  basing  whatever  further  analysis  is  desired  on  that 
model’s  established  properties;  this  aspect  is  more  fully  discussed  In 
Section  2.  A  further  point  is  that  our  focus  is  on  current  and  projected 
trend  rather  than  on  extracting  patterns  in  historical  series.  While  ex 
post  seasonal  adjustment  and  trend  estimation  are  important  for  some  pur¬ 
poses,  the  problem  of  prime  importance  to  many  people,  including  forecasters, 
planners  and  policymakers,  is  obtaining  the  best  estimate  of  what  is  happen¬ 
ing  now.  Finally,  the  procedure  is  determined  uniquely  by  the  model  and 
thus  avoids  the  elements  of  arbitrariness  inherent  in  stochastic  modelling 
approaches  such  as  those  of  Box,  Hillmer,  and  TIao  (1978)  and  Pierce  (1978). 

Section  2  of  the  paper  discusses  the  importance  of  model-building 
in  scientific  inquiry  and  introduces  the  time  series  models  on  which  the 
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present  investigation  is  based.  In  Section  3  these  models  are  used  to 
develop  estimates  of  trend  for  general  seasonal  series.  This  involves  an 
analysis  of  the  forecast  function  of  the  modelled  series,  expressed  as  the 
solution  of  a  difference  equation.  This  solution  contains  a  set  of  adaptive 
coefficients  (adaptive  to  the  current  time  period),  some  of  which  can  be 
associated  with  the  tries'  current  trend  (or  growth  rate  for  logged  series). 

Section  4  studies  modelling  and  trend  estimation  for  the  series 
of  demand  deposits  at  commercial  banks,  a  major  component  of  the  U.S.  money 
supolv.  Estimated  and  actual  growth  rates  are  compared  over  a  5-year  period, 
and  their  mean  square  error  is  found  to  correspond  closely  to  the  demand- 
deposit  model's  error  estimate.  Some  concluding  remarks  comprise  Section  5. 
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2.  TIME  SERIES  MODELS  AND  MODELLING  PHILOSOPHY 


Seasonal  adjustment*  trend  estimation  and  so  forth  imply  a  belief 
in  regularities  of  some  kind  in  the  series  studied.  This  "sequence  regular¬ 
ity"  is  imperfect  and  therefore  must  be  expressed  in  terms  of  probability. 
Sequence  regularity  expressed  by  probability  is  precisely  what  a  time  series 
model  is.  While  a  time  series  model  implies  a  particular  choice  of  procedure, 
a  particular  choice  of  procedure  implies  a  particular  kind  of  time  series 
model.  Thus,  it  was  shown  by  Cleveland  and  Tiao  (1976)  that  the  Xll  program 
implies  (very  nearly)  a  particular  kind  of  time  series  model  with  particular 
values  of  its  adjustable  parameters.  Let  us  call  this  model  M(X11).  Since 
Xll  was  arrived  at  by  highly  skilled  people  using  many  iterations,  M(X11)  is 
an  average  or  compromise  model  for  the  range  of  series  on  which  its  organizors 
tested  it.  Consequently,  Xll  can  be  expected  to  do  an  averagely  good  job  for 
series  whose  models  are  similar  to  M(X11)  (for  series  which  are  like  the  ones 
it  was  tested  on)  and  not  as  good  job  on  other  series.  For  example,  if  Xll 
were  applied  to  a  random  series,  it  would  induce  seasonality  in  the  series. 

The  above  implies,  of  course,  that  models  and  methods  are  not  arbi¬ 
trary  but  can  be  built  and  applied  in  a  series  of  logical  steps.  These 
steps  have  been  set  out,  for  example,  by  Box  and  Jenkins  (1970).  Basically, 
the  argument  is  that  a  model  is  a  transformation  of  data  to  white  noise, 
evidenced  by  residuals,  uncorrelated  with  any  other  known  input.  The  model 
building  process  is  an  iteration  which  is  guided  at  each  step  by  the  need  to 
achieve  this. 

The  conclusion  from  this  argument  is  as  follows  *  Since  different 
time  series  have  different  probability  structures  (i.e.,  different  models), 
methods  of  analysis  such  as  seasonal  adjustment  methods,  trend  estimation 
methods,  etc.,  should  be  different  and  should  depend  on  the  model.  Similarly, 


seasonal  adjustment  methods,  if  they  are  really  relevant  to  the  problem,  must 
be  model-based  [as  are  those  proposed  by  Box,  Hillmer,  and  Tiao  (1978)  and 
Pierce  (1978),  and  Burman  (1980)],  so  that  they  adapt  to  the  series  under  study. 

In  any  statistical  or  scientific  investigation,  therefore,  one 1 s 
approach  should  consist  of  two  stages: 

1)  to  build  a  model  for  the  data  under  study; 

2)  to  use  that  model  to  supply  answers  to  whatever  it  is  that 
we  want  to  know. 

This  process  is  necessarily  ongoing,  as  seldom  is  a  model  ever  ideal;  rather, 
model  building,  like  the  pursuit  of  happiness,  is  something  that  one  should 
always  be  working  on,  continually  examining  residuals  and  worrying  about 
other  variables. 

In  the  present  study,  we  construct  seasonal  ARIMA  models  for  time 
series  and  show  how  they  can  be  used  to  estimate  the  trend  (or  growth  rate) 
of  the  series.  While  the  various  aspects  of  a  series1  behavior  quite  likely 
result  from  forces  not  entirely  captured  by  ARIMA  models,  we  believe  that 
thi9  is  a  natural  place  to  begin.  Having  first  constructed  an  ARIMA  model, 
we  can  then  allow  for  outliers  or  interventions,  or  relate  two  or  more  time 
series.  The  need  for  this  more  elaborate  analysis  is  suggested  both  by  an 
examination  of  the  white  noise  residuals  —  the  sequence  of  one-step- ahead 
forecast  errors  —  and  by  economic  or  institutional  knowledge  of  any  addi¬ 
tional  variables  to  which  the  residuals  may  be  related.  The  criterion  for 
Including  the  results  of  such  an  analysis  is  whether  the  variance  of  these 
residuals  is  materially  reduced. 

An  ARIMA  model  is  a  means  for  expressing  the  current  value  zt  of  a 
time  series  in  terms  of  the  past  values  zt_i ,  Zt-2 *  •  •  •  that  series.  Thus 


we  could  write 


(2.1) 
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Zt  -  n*t-l  +  *2zt-2  +  ...  +  ac  -  2>jZt-j  +  at 

where  at  is  the  random  error  (shock,  innovation)  at  time  t  unpredictable  from 
the  series*  past,  assumed  to  be  (iid)  with  mean  0  and  variance  .  If  B 
denotes  the  backshift  operator  defined  by  B^zc  *  z^-j  then  (2*1)  becomes 

(1  -  »  tt(B)zc  »  at  (2.2) 

As  this  expression  implies  a  potentially  infinite  number  of  parameters 
we  assume,  for  parsimony  and  simplicity,  that  tt(B)  can  be  expressed  as 
a  ratio  <&(B)/9(B),  whereupon  multiplying  (2.2)  through  by  6(B)  we  obtain 

$(B)zt  *  <KB)at.  (2.3) 

In  practice  zt  is  frequently  nonstationary  but  stationary  after  differencing 
the  series  d  times,  in  which  case 


♦(B)  -  A(B)  0(B)  (2.4) 

where  A(B)  »  (l-B)^.  If  <j>( B)  and  9(B)  are  of  degrees  p  and  q  then  (2.3)  is 
referred  to  as  an  autoregressive-integrated  moving  average  (ARIMA)  model  of 
order  (p,d,q).  In  the  model  (2.3),  the  roots  of  the  characteristic  equations 
4>( B )  *  0  and  9(B)  •  0  lie  outside  the  unit  circle,  and  the  zeroes  of  the 
Mdif f erencing"  polynomial  A(B)  lie  on  the  unit  circle.  The  differenced 
series  wt  *  A(B)zt  is  therefore  a  stationary,  invertible  autoregressive- 
moving  average  processes 

p  q 

wt  »  Z  <f>4Wt_4  -  Z  04at-4  4*  at.  (2.5) 

j-1  J-l 

In  constructing  a  model  for  an  observed  time  series  it  is  essential 


to  have  available  means  for  identification  (of  the  particular  model  form) ,  for 
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estimation  of  this  form,  and  for  diagnostic  checking,  or  model  criticism, 
after  fitting.  For  time  series  this  procedure  is  discussed  in  Box  and  Jenkins 
(1970,  Ch.  6-8). 


* 

* 


3.  SEASONAL  MODEL  AND  TREND  ESTIMATION 


As  indicated  in  Section  2,  a  model  for  a  series  or  system  which 
captures  the  essential  features  of  that  process  may  be  used  to  provide  answers 
to  the  major  questions  of  this  investigation.  Here  we  are  interested  in  the 
use  of  a  model  fitted  to  a  time  series  to  construct  estimates  of  that  series’ 
current  trend  or  rate  of  growth.  These  quantities  will  be  seen  to  depend  on 
forecasts  of  the  series’  future  values,  and  thus  this  section  begins  by  ana¬ 
lyzing  forecast  functions  for  ARIMA  models,  particular  for  seasonal  ARIMA 
models,  extending  Chapters  5  and  9  of  (Box  and  Jenkins,  1970).  We  then  pre¬ 
sent,  in  Section  3.3,  the  general  trend  estimation  procedure. 

3 .1  Forecast  Function  for  ARIMA  Models 

Given  a  segment  zt,  zt_i , . . .  of  a  time  series,  the  forecast  of  a 
future  value  of  that  series  with  minimum  mean  square  error  is  the  con¬ 

ditional  expectation 

z  £  (  £ )  *  Et  ( z  m  E  ( z  £  I  z  £  t  (3.1) 

For  the  model  (2.3)  this  forecast  function  (considered  as  a  function  of  £, 
the  '‘lead  time”)  is  known  to  satisfy  the  difference  equation 

$(B) zt ( £)  -  0(B)at(£)  (3.2) 


where  B  now  operates  on  l,  and 


In  particular, 


a  t(£) 


at+£  *  A  <  0 

0,  l  >  0. 


<&( B) zt( £)  »  0 ,  1  >  q . 


(3.3) 
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For  l  >  q  -  P,  eq.  (3.2)  has  a  general  solution  of  the  form 

P 

Zt(*)  -  Z  (3.4) 

i-1 

where  the  f £ ( A >  *  1  _<  i  £  P,  are  functions  (in  general  polynomials,  exponen¬ 
tials,  sines  and  cosines  and  products  of  these)  of  the  lead  time  l;  in 
particular,  f^(£)  »  g^  if  the  root  of  $(B)  *  0  is  real  and  distinct. 

Having  determined  the  functions  f  £  ( Jt)  ,  the  solution  (3.4)  further 
requires  a  knowledge  of  the  adaptive  coefficients  so  named 

because  they  change  with  the  forecast  origin  t,  adapting  to  new  information 
(data)  continually  becoming  available.  In  Appendix  1  it  is  shown,  in  turn, 
how  to  obtain  the  {b£c^}  : 

(a)  from  an  Initial  set  of  forecasts, 

(b)  directly  from  the  available  observations  {zt, 

(c)  from  their  previous  values  as  additional  observations 
become  available  and  the  forecast  origin  is  shifted 
forward . 

3 . 2  Seasonal  ARIMA  Models 

An  important  case  of  the  ARIMA  model  (2.3)  found  useful  In  analyzing 
seasonal  time  series,  which  exhibit  periodic  behavior  of  a  stochastic  or  adap¬ 
tive  nature,  is  the  multiplicative  seasonal  model.  Suppose  that  in  (2.3)  the 
operators  $(B)  and  8(B)  factor  according  to 

*(B)  -  «1CB)«2(Bs) 


3(B)  -  el(B)B2(Bs) 


where  the  seasonal  operators  $2  and  ^  are  functions  of  Bs,  s  denoting  the 
period  (e.g.,  s*12  for  monthly  data  with  an  annual  seasonal  pattern).  Then 
the  model  (2.3)  may  be  written 

<t1(B)4>2(Bs)zc  =  91(B)02(Bs)at.  (3.5) 

A  simple  example  of  a  seasonal  ARIMA  model  would  be 


zt  “  zt-12  +  at  ^t-l  » 
or 

(1  -  B12)zt  -  (1  -  0B)at, 

in  which  ^(B)  -  1  -  B12  ,  9(B)  -  1  -  9B ,  and  ^(B)  -  92(B)  *  1.  This  model 
expresses  the  current  observation  as  the  sum  of  last  year's  value  and  a  linear 
combination  of  this  and  last  month's  shocks. 

Typically  and  9^  will  be  of  low  order.  As  in  (2.3),  ^(B) 
incorporates  nonseasonal  differencing  operators  applied  to  render  the  series 
stationary.  For  example,  the  model  fitted  to  the  Demand  Deposit  series  in 
Section  4  is  of  the  form 


(1  -  B) ( 1  -  BI2)zt  -  (1  -  9B) ( 1  -  0B12)  (3.6) 

and  has  [noting  1  -  B^-2  =*  (1  -  B)(l  +  B  4*  ...  +  B^)] 

11  . 

♦L(B)  -  (1  -  B)2,  *2<B)  -  ?.  BJ.  (3.7) 

j=0 


3 .3  Estimated  Trend 

The  general  procedure  for  trend  estimation  for  a  seasonal  time 
series  can  now  be  set  forth.  In  the  seasonal  ARIMA  model  (3.5),  if  ,  $2 , 

9|  and  09  are  respectively  of  degrees  P^,  P2 ,  qi  and  q2  B,  then  analogous 
to  equation  (3.4)  the  eventual  forecast  function  for  the  model  (3.5)  may  be 


written  in  the  form 


in 


Zt(i)  -  r  b{jr>  fua)  +  E  b^>  f21(l) 

1-1  1-1 

-  Tt(£)  +  St(£)  (3.8) 

The  second  sum  on  the  right  hand  side  of  (3*8)  is  associated  with 
the  seasonal  par*-  of  the  model,  determined  by  #2  an<*  •  The  f*rsC  sum  Tt(0 
embodies  all  systematic  nonseasonal  influences,  and  may  thus  be  regarded  as 
the  trend-  Our  proposed  trend  estimate  for  seasonal  time  series  described 
by  the  model  (3*5)  is  therefore 

P 

1 

Tt(t)  -  *t<£)  -  St(£)  -  E  b<p  fu(t)  (3.9) 

i-1 

The  quantity  Tt(£),  which  incorporates  anything  systematic  and 
nonseasonal  that  is  known  about  the  series  at  the  current  time  t,  will  often 
correspond  well  to  our  instinctive  notion  of  trend.  For  example  if  $(B) 
contains  a  factor  (1  -  B)^  reflecting  a  nonstationari ty  removed  by  differenc¬ 
ing  the  series  d  times,  then  Tt(£)  contains  a  polynomial  of  degree  d-1  (in 
i) ,  with  coefficients  adapting  with  the  origin  t.  Thus  if  d  is  greater 

than  1  (equal  to  or  greater  than  1  if  a  constant  term  also  appears  in  the 
expression  for  (1  -  B)zt)  then  at  least  a  linear  trend  will  be  incorporated 
into  Tt(£).  On  the  other  hand,  we  would  not  ordinarily  think  of  a  series 
zt  which  is  already  stationary  as  possessing  a  trend,  and  indeed  Tt(£)  for 
such  a  series  would  either  be  zero  or  possess  only  damped  terms  based  on  the 
roots  of  the  stationary  autoregressive  operator  40). 

Similar  statements  can  be  made  about  St(Jl);  for  example,  whenever 
is  a  factor  of  St(£)  will  include  a  periodic  function  of  period 

s;  this  is  the  case  for  the  model  (3.6). 


-  XI  - 

4.  ANALYSIS  OF  DEMAND  DEPOSIT  SERIES 
The  U.S.  Money  Supply  is  a  closely  watched  and  much  investigated 
series  in  economic  policy  making  and  analysis.  Changes  in  the  underlying 
behavior  of  this  series  have  an  immediate  impact  on  financial  markets  and 
are  felt  to  have  a  strong  association  with  many  measures  of  economic  activity, 
including  income,  inflation  and  employment.  It  is  therefore  of  prime  impor¬ 
tance  for  both  the  Federal  Reserve  Board  and  the  public  to  have  available  the 
best  and  most  timely  information  possible  on  the  current  behavior  of  this 
series,  i.e.,  its  current  trend  or  rate  of  growth.  A  key  element  of  the 
money  supply  is  its  demand-deposit  component,  and  this  section  is  devoted  to 
estimating  the  current  rate  of  growth  of  this  series. 

4 . 1  Airline  Model 

The  model  fitted  to  the  Demand  Deposit  series  (Section  4.2)  is  of 
the  general  form  (3.6),  which  has  been  found  to  describe  quite  well  a  number 
of  seasonal  time  series.  It  was  first  fitted  by  Box  and  Jenkins  (Chapter  9)  to 
a  series  of  logged  monthly  passenger  totals  in  international  air  travel  and 
has  thus  become  known  as  the  Airline  model.  This  model  has  a  forecast  func- 
tio  zt(£)  satisfying  the  difference  equation 

(1  -  B)(l  -  B12)zt(*)  -  (1  -  9B) ( 1  -  0B12)at(O  (4.1) 

where  at(£)  is  as  below  (3.2).  This  equation  has  a  solution  which  can  be 
expressed  in  several  equivalent  forms,  involving  dummy  variables  or  sines  and 
cosines  for  the  seasonal  term  St(£)  in  (3.8).  These  are  derived  in  part  by 
Box  and  Jenkins  (1970)  and  in  part  in  Appendix  2,  where  it  is  found  that  the 
most  appropriate  expression  for  our  current  use  is  of  the  form 


zt(0  -  +  b^  +  b(c)*  . 


(4.2) 
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The  three  terms  on  the  right  hand  side  of  (4*2)  are  respectively  the  level, 
seasonal  and  slope  components  of  the  forecast.  Therefore  the  estimated 
trend  is,  as  in  (3.9), 

Tt(l)  -  zt(*)  -  b<c)  =  b<c)  +  b(c)*  (4.3) 

where  in  (4.2)  and  (4.3)  m  »  1,2 ,..,12  indexes  the  month  of  the  year  in  which 
the  value  zt+£  falls.  There  are  thus  14  coefficients  (b£ , . . .  ^  > 

b^c^)  the  first  12  of  which  satisfy  the  constraint  £  b£c^  31  0 .  As  with  the 
coefficients  in  nonseasonal  ARIMA  models  these  may  also  be  determined  In  several 
ways,  as  illustrated  in  Appendix  2.  In  particular  the  following  results  are 
of  primary  interest: 

1.  The  slope  coefficient  in  (4.2)  is  given  by 

b(t)  -  [zt(13)  -  zc(l)]/12,  (4.4) 

which  is,  as  perhaps  expected,  the  average  monthly  increment  to  the  forecasts. 
The  other  b-coef f Icients  are  similarly  determined  from  the  forecasts  zt(l ),..., 
zt(13) . 

2.  The  b’s  are  linear  combinations  of  forecasts  which  are  them¬ 
selves  linear  combinations  of  the  data  zt,  Zfl » •  •  •  •  Thus  the  coefficients 
b,  bQ,  , . . . ,  b]^  are  expressable  directly  in  terms  of  this  available  data, 
as 

b[c)  =■  W1(B)zt.  (4.5) 

For  example,  Figure  1  shows  the  weights  |w j  }  in  determining  the  slope 

b(t)  -  W(B)zc  -  IW,zt_.,  (4.6) 


for  the  Airline  data  referred  to  above. 
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3.  The  coefficients  of  the  forecast  function  evolve  as  new  obser¬ 
vations  result  in  a  shifting  forward  of  the  forecast  origin,  t;  for  example, 
the  relation 

b(c+1)  =  b(c)  +  (1  -  9)0  -  0)at+1  (4.7) 

shows  how  a  fraction  of  the  new  information  (innovation,  shock)  at+i  is  incor¬ 
porated  into  the  revised  slope  estimate, 

4 . 2  ARIMA  Model  for  Demand  Deposit  Series 

Monthly  observations  on  the  demand-deposit  component  of  the  money 
supply  (DD)  were  obtained  for  the  period  1/69  -  1/78  (Source:  Federal  Reserve 
Bulletin,  The  data  are  periodically  revised,  the  series  here  used  being  as 
of  10/78).  Application  of  the  usual  model  fitting  procedure  (Box  and  Jenkins. 
1970,  Ch.  6-8)  to  the  logged  series  ( 3  log  DDt)  yielded  a  fitted  ARIMA 
model  of  the  form 

AA17zt  »  (1  -  -6B12 )( 1  +  .lB)at  (4.8) 

with  aa  31  .0058  .  This  is  the  Airline  model  form,  with  parameter  estimates 
9  »  .6,  9  *  -.1,  so  that  the  trend  in  the  logged  demand  deposit  series  can  be 
estimated  as  in  Section  3.3.  We  shall  focus  on  the  first  difference  of  this 
trend,  which  is  the  estimated  rate  of  growth  ( deseasonalized  or  seasonally 
adjusted)  of  demand  deposits.  It  is  In  fact  the  rates  of  growth  of  money 
supply  series  which  are  most  often  examined  both  within  the  Federal  Reserve 
and  in  financial  communities. 

4 .3  Estimated  Rate  of  Growth 

Taking  the  first  difference  of  the  trend  function  (4.3)  shows  that 
the  estimated  current  (at  time  origin  t)  rate  of  growth  of  demand  deposits  is 
simply  the  slope  coefficient  b(c).  Thus,  by  virtue  of  eq .  (4.4)  the  practical 
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application  of  the  concepts  herin  proposed  reduces  in  this  case  simply  to 
computing  the  difference  between  the  lead-13  and  lead-1  forecasts  of  log  DD. 

As  it  is  customary  to  express  these  figures  in  terms  of  annual  rates,  we 
omit  the  division  by  12. 

Figure  2  shows  the  weights  which  can  be  used  to  calculate  b^c^ 
from  the  series  values  zt,  as  in  equation  (4.6)  (without  the  divi¬ 

sion  by  12  in  (A2.ll)).  This  is  broadly  similar  to  Figure  1  for  the  airline 
data,  differences  resulting  from  che  first  order  MA  parameter. 

The  calculation  of  the  estimated  growth  rate  for  DD  was  made 

over  the  period  8/72  -  9/78,  and  the  results  are  shown  in  Table  1  and  Figure  3. 

4,4  Error  Limits  of  Growth-Rate  Estimates 

Since  the  slope  estimate  and  the  other  coefficients  b^^  are 

functions  of  the  forecasts,  the  errors  in  these  adaptive  coefficients  are  the 
corresponding  functions  of  the  forecast  errors  et(£)  *  Zt+£  “  zt(£)«  In  par¬ 
ticular  the  error  corresponding  to  b^c^  in  (4.4)  is 

e(t)  -  [et(13)  -  et(l)] 

“  t at+13  +  ^  at+12  +  +  'I'll  at+2 

+  (i>  -  1)  at+1]  (4.9) 

12 

with  mean  square 

11 

a2  [eft>  ]  -  [£0  4^  +  ( ’h.2  *  D2]  •  (4.10) 

Letting  the  actual  growth  rates  be  defined  by 

8(0  *  zt+13  ~  zc+l 

also  gives 


e(  t)  rn  b(  c )  • 
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Table  1  and  Figure  3  display  these  actual  growth  rates  and 

errors  The  dotted  lines  in  Figure  3  are  standard  error  limits  ±  e(c) 

associated  with  the  estimated  rates  of  change.  Based  on  (4.10)  we  have 

E[bC C)  -  0(t)j2  »  .00040,  (4.11) 

making  the  standard  error  .020  or  2%.  (Recall  that  these  are  annualized 
figures;  the  monthly  standard  error  itself  is  (1/12 )th  of  this  amount.)  This 
mean  square  error  was  obtained  from  the  fitted  model  (4.8),  and  indeed  one  of 
the  virtues  of  the  model  building  approach  is  that  such  internal  error  limits 
are  thereby  available.  But  it  is  also  of  interest  that  the  empirical  MSE 
over  this  period, 

-  g(t)]2  =,  .00044,  (4.12) 

6 1 

is  in  close  agreement  with  the  model-based  estimate. 

Figure  3  also  illustrates  the  limitations  of  the  present  model  and 
the  importance  of  making  the  model  building  effort  an  ongoing  one,  to  incor¬ 
porate  additional  relevant  information.  To  give  an  estimate  of  (say)  6%  as 
the  current  seasonally  adjusted  rate  of  increase  of  demand  deposits  but  to 
have  the  95%  confidence  interval  span  the  range  2  to  10%  is  uncomfortable, 
even  though  other  money  supply  estimates/ forecasts  typically  have  errors  of 
comparable  magnitudes  [e.g.,  Porter  et  al.  (1978),  Pierce  et.  al.  (1981)]. 

Ths  use  of  other  variables  related  to  causes  of  the  shortfall  of  money  growth 
during  1974  (relative  to  its  estimate  in  Figure  3)  or  of  the  overshooting  of 
money  growth  in  1976-77  could  possibly  reduce  this  error. 
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5.  CONCLUSIONS 

We  have  presented  a  method  for  estimating  the  current  trend  or  rate 
of  growth  of  seasonal  time  series,  based  on  models  for  those  series,  and  have 
illustrated  its  application  with  the  demand  deposit  component  of  the  U.S. 
Money  Supply.  The  resulting  trend  estimates  are  thus  derived  in  an  optimal 
manner  from  the  underlying  characteristics  of  the  series  and  consequently  do 
not  suffer  from  the  arbitrariness  inherent  in  methods  not  so  derived.  But  we 
have  also  noted  the  need  to  regard  such  model  building  investigations  as 
ongoing  pursuits;  in  particular,  the  incorporation  of  outliers,  interventions 
and  other  variables  into  the  model  should  reduce  the  mean  square  error  of 
the  resulting  trend  or  growth-rate  estimates. 
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APPENDIX  1.  ADAPTIVE  COEFFICIENTS  IN  ARIMA  FORECASTS 
AND  TREND  ESTIMATES 

A. I  Determination  of  Adaptive  Coefficients  from  Initial  Forecasts 

For  lead  times  Z  >  q,  the  forecast  function  2^(0  satisfies  the 
homogenous  difference  equation  (3.3),  and  thus  for  such  Z  (3.4)  holds  for 
arbitrary  constants  These  constants  are  determined  (uniquely)  by  a 

set  of  initial  conditions  which  are  equivalent  to  assuming  that  (3.4)  [though 
not  (3.3)]  also  holds  for  Z  =*  q-P+1 ,  ...,  q-1 ,  q,  which,  given  zt(  Z)  for 
these  Z  values  [z t(Z)  =*  zt+£  If  Z  £  0]  ,  defines  a  set  of  P  equations  in  the 
P  unknowns  b£c\  b^^.  In  fact,  any  set  of  P  consecutive  forecasts/ 

actuals 

z (Z)  «  [z t(Z),  zt(£)+l),  ...,  zt(*+P-l)r  (Al.l) 


such  that  Z  >  q-P  can  be  used.  Letting 


b(t)  ■  [b(c^  ,  ]  and 


fl  U) 


fx(  JW-P-1) 


fp(  O 


fP( t+p-l) 


( A1 . 2 ) 


this  system  can  be  written  In  the  form 


zt(l)  =■  Ftb(c), 


whence 


b<c> 


f  A1 .  3  ) 


(A1.4) 
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If  q  _<  P ,  it  is  often  convenient  to  consider  the  system  in  terms  of  che  first 
P  forecasts  (£  38  1),  so  that,  with  88  z^, 


k(  t )  p  1  1 

Z  1,1  If 


(A1.5) 


As  an  example,  consider  the  AR(1)  process 


for  which  (3.4)  is 


so  that 


(1  -  $B)xt  3  at 


Zc(t)  =  f (l)  b(c)  =  <t>lzt,  l  >  0, 


and  the  calculation  (Al.4), 


F,  -  f(l)  -  , 


*  (6^)~*  zt(£)  3  zt» 


can  be  made  for  any  l  >  q-P  =  *1 , 


Al.2  Adaptive  Coefficients  in  Terms  of  the  Observed  Series 

Since  the  coefficients  are  linear  combinations  of  an  initial 

set  of  forecasts  [Eq.  (Al.4)],  and  since  the  forecasts  at  origin  t  are  linear 
combinations  of  the  series  (zt,  , . .  •  \  (Whittle,  l%3,  Box  and  Jenkins, 

L970),  it  follows  that  the  coefficients  are  expressable  as  linear  combinations 
of  the  series.  In  particular,  writing 


zt(  l) 


(A1.6) 


It.  follows  chat 


ir'  *■'  ( B) 


"(B)  [  Z  tp.BJ  el, 
.1- 1 


(Al.7) 
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where 


Hence,  if 


it(B)  -  r  TTjBJ  =»  <jT 1  ( B )  =«  (  Z  'i'jBJ)"1 


tt(^)(B)  -  [ir(*)(B) .  ff(Z+P-l)(B)]  '  (A1.8) 


then  from  (A1.4)  and  (A1.6),  the  desired  expression  for  the  coefficients  is 


-  W(B)zc 


(A1.9) 


where 

W(B)  -  F^1  TT^^CB)  -  [W1(B),...,Wp(B)l  '  •  (A1.10) 

Of  course  the  weights  W(B)  are  the  same  for  any  value  of  Z  >  q-P  in  (Al.10). 
For  the  AR(1)  process  introduced  at  the  end  of  Section  3.1, 


tt<  (B>  =  <t>z,  W(B)  =  1. 


A1 .3  Updating 

Given  a  new  observation  zt+^,  the  forecasts  from  origin  t+1  may  be 
expressed  as  simple  adjustments  to  the  forecasts  calculated  at  origin  t,  the 
adjustments  depending  on  (and  only  on)  the  new  information  at+i*  Thus  the 
adaptive  coefficients  may  be  so  updated.  In  particular,  the  usual  updating 
formula 

zt+l(*>  -  Zt(m)  +  «t+l 

Is,  In  the  previous  notation. 


;(£+l)  +  ,i.  a 

t  h  at+r 


(Al.ll) 


with  i i  •  •  • » ^£+p-i )  >  so  that,  from  (A1.3)  and  its  counterpart 


at  t ime  t+1 , 
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whence 


It  b(C+U  =  l2+lbU)  +  ^ac+l 

u(  C+l  )  a  p“l  -p  w(t)  +  t?-!  -Ji  a 

Z  I*  ZMz  F*  Vt+i 

=  L  b^^  +  k  a. 


(A1.12) 


t+I* 

Again  considering  the  AR(1)  process  at  the  end  of  Section  3.1, 
(Al.13)  becomes,  taking  Z  *  0  so  that  FQ  =»  1  and  F|  *  <p, 


(A1.13) 


,(t+n  =  *h(t) 


<(>b^L'  +  a 


t+1 


which,  since  =  zt,  is  simply  Che  equation  for  the  process. 

A1 .4  "Error*’  of  Adaptive  Coefficient  Estimates 

Since  the  various  b(t)fs  are  functions  of  the  forecasts,  the  errors 
in  these  adaptive  coefficients  may  be  defined  as  the  corresponding  functions 
of  the  forecast  errors.  In  general, 


e[b^ch  =  F_1  e£^  -  _  h(c) 


CA1 .14) 


where,  analogous  to  (Al.4), 


=*  FjJ  (zt+i ,  •  •  * ,  2t+t+p-i) 


(Al.15) 


Note  that  and  efb^)],  unlike  b(^)  itself,  depend  on  the  value  of  Z. 


1 

1 
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APPENDIX  2:  AIRLINE  MODEL  TREND  ESTIMATION 


The  9olution  of  the  forecast  function  (4.1)  for  the  airline  model 
(3.6)  is  of  the  form 


zta> 


2  ttJ  5, 

sin  - ]  . 

12 


An  equivalent  form,  which  is  easier  to  work  with,  uses  monthly  dummy  or  indi¬ 
cator  variables,  and  there  are  several  variants.  Adapting  the  development  in 
[Box  and  Jenkins  (1970),  p.  309],  it  can  be  shown  that  (A2.1)  can  be  written 
in  the  alternative  form 


zt(i) 


b^;)  +  b(t)£ 
om 


(A2.2) 


where  l  is  measured  in  months. 

Next,  note  that  the  monthly  coefficient  incorporates  both 

an  overall  level  effect 


»it> 


_1 

12 


12 

l 

m-1 


b(t) 

om 


(A2.3) 


and  a  seasonal  effect  specific  to  the  m month 


,(t)  ,  v(t)  .  v(t), 


b^U'  -  b^ 
om  o 


(A2.4) 


thus,  the  forecast  (A2.2)  is  equivalently  equation  (4.2)  of  the  text,  from 
which  the  trend  (4.3)  is  derived. 

Parallel  to  Appendix  1  we  examine  each  group  of  coefficients,  obtain¬ 
ing  expressions  for  them  (a)  from  an  initial  set  of  forecasts,  (b)  from  the 
available  observations  zt,  and  (c)  from  their  previous  values  as  new 

observations  become  available  and  the  origin  is  shifted  forward. 
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A2 . 1  Initial  Determination  of  Adaptive  Coefficients 

As  in  (Al.l)  let  denote  the  vector  of  P*13  forecasts  zt(£), 

■  zt(£+12),  and  for  convenience  let  us  take  t  *  1 ,  letting  *  zt, 

though  any  positive  l  would  suffice.  Then  (A1.3)  is  a  system  of  13  equations 

of  the  form  (A2.2)  which  needs  to  be  solved  for  the  13  unknowns  »a> . 

bo^2'  b<°-  Thls  is  most  immediately  done  by  noting  that  the  equations 
for  zt(l)  and  zt(13)  imply  that  the  slope  coefficient  is  equation  (4.4)  of 
the  text. 

In  terms  of  the  other  12  equations  then  each  have  solutions 

bom^  ”  ztO)  ”  m  b^  (A2.5) 

giving  the  mc^  monthly  coefficient  as  an  adjustment  to  the  lead-m  forecast, 
the  correction  being  the  amount  required  to  offset  the  trend  component  of 
the  forecast. 

Given  (4.4)  and  (A2.5),  the  level  and  seasonal  effect  coefficients 
b^c^  and  b^c^  in  (4.2)  are  determined  as  in  equations  (A2.3)  and  (A2.4): 

12 

b<c)  =  _L  E  z  (m)  -  11  [z't(13)  -  zt(l)],  (A2 .6) 

12  m*!  24 

12 

b£c)  -  it(«)  -  _L  Z  zt(n)  +  lljli™  [£t(l3)  -  zt(ni  •  (A2.7) 

12  n-1  C  24  c  c 

A2 .2  Adaptive  Coefficients  in  Terms  of  Available  Data 

As  In  Appendix  1,  the  various  adaptive  coefficients  are  each  expres- 
sable  as  linear  combinations  of  the  observed  series  zc,  zc_^,...  .  Considering 
first  the  forecasts  written  as  In  (A2.2)  containing  12  monthly  levels  , 
let 

b^  -  W0m(B>zc,  12  (A2.8) 


t 
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Then  from  (Al.6)  and  (A2.5), 

Wom(B)  ”  (B>  -  ir(l)(B)];  (A2.9) 

and  writing  the  slope  coefficient  as 

b(t)  =  W(B)zt  (A2.10) 

which  is  eq.  (4.6)  of  the  text,  it  follows  that 


W(B)  -  _L  [ttU3)(B)  _  tt(1)(b)].  (A2.ll) 

12 

Let  us  next  express  the  forecast  function  as  (4.2),  in  which  the  monthly 

level  is  separated  into  an  overall  level  b^^  and  a  seasonal  effect 

o  m 

(1^1  b<c>  -  0).  From  (A2.6)  and  (A2.7),  If 


b,(£> 


-  VB)zt»  ID"0-  1’"*> 12 


(A2.12) 


then 


W0(B)  -  -i  E  Wom(B) 


12 


_L  E  ir(®)(B)  -II  [  ttC13 ) (B)  -  *<1)(B)] 
12  24 


(A2.13) 


and 


Wm(B)  -  Wom(B)  -  W0(B),  m-l,...,12 


12 


TT(m)(B)  -  _1  z  tt(*0(b)  +  13  ~  2m  (ir(13)(B)  -  *<1)(B)1. 
12  n-l  24 


(A2.14) 


Whereas  the  transformation  from  the  forecasts  zt  to  the  adaptive 
coefficients  b(0  does  not  depend  on  the  airline-model  parameters  (as  the 
form  of  eventual  forecast  function  does  not),  the  forecasts  themselves, 
and  thus  the  weights  such  as  W(B)  and  Wm(B),  are  functions  of  9  and  0. 


28  - 


These  weights  functions  can  be  determined  for  the  values  9  *  .4  and  0  *  .6, 
corresponding  to  the  least  squares  fit  obtained  by  BJ  for  the  airline  data. 

Using  numerical  recursive  calculations v  we  first  determine 

ir(B)  =  ckB) 

-  (1-B)  (1-B12)  (1-.4B)-1  (1-.6B12)"1  CA2.15) 

truncated  at  order  256,  and  then  ff(0(B)  for  l  M  l ,...,13  is  obtained  as  in 
(Al.7).  Next  we  form  the  various  weight  functions,  including  W(B)  which  was 
plotted  in  Figure  1.  Note  the  periodicity  exhibited  by  this  weight  function, 
which  was  true  of  the  others  as  well. 


A2 .3  Updating 

As  discussed  in  A1.3,  the  coefficients  of  the  forecast  function 
(3.4)  evolve  as  new  observations  result  in  a  shifting  forward  of  the  forecast 
origin.  The  basic  formula  for  the  calculation  is  equation  (Al.13),  with  P=*13. 
However,  it  is  easier  to  work  directly  with  the  airline  model  form,  adapting 
the  development  in  [Box  and  Jenkins  (1970),  pp.  310-11]  to  the  modified  form 
(A2.2)  of  the  forecast  function. 

Letting 


1  =  1-9,  A  -  1  -  0, 
and  letting  1  31  12r  +  m,  1  _<  m  <_  12,  then 


'H  *  >h2r+m  *  *(!  +  rA)  +  512m 


(A2.16) 


where  $\2m  *  1  Iff  m  *  12  and  0  otherwise.  From  (A2.16)  and  the  updating 
formula  (Al.13),  it  can  be  shown  that 

b(t+1)  -  b<c>  +  (AA/12)at+1, 

bU+1)  ■  *  b<t)  *  ^  +  su.‘l*t+f 


(A2.17) 


(A2.1S) 
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where,  if  m  ■  12,  is  replaced  by  * 

Thus  the  previously  forecasted  series  level  for  each  month  is 
adjusted  by  the  previously  forecasted  increment  to  the  trend,  now  incorporated 
into  the  level  at  the  new  origin;  and  all  the  coefficients  are  adjusted-  by 
varying  fractions  of  the  new  information  (innovation,  shock)  a^+1 * 

For  the  airline  model  coefficients  X  *  .6  and  A  -  .4,  (A2.17)  and 
(A2.18)  become 


b^t+1^  *  +  .02 


‘t+l 


bom+1)  *  bnCl+ 1  +  b(C)  + 
om  o ,  urri 


(.6  -  .02m)  at+*j_ ,  m  ^  12 
.76  at+i ,  m  -  12 


The  seasonal  effects  and  the  level  effect  in  (4.2)  are 

in  o 

updated  in  corresponding  ways,  based  on  (A2.3)  and  (A2.4).  Thus, 


^(t+l)  »  b^^  +  +  [X(l  -  12  A)  +  A]  at+^> 


24 


(A2.19) 


A 


since  (Em)/144  -  13/24  and  and 

bitM)  -  +  [  -  <l  -  512»>A1  *t+l  <«-20> 


by  subtracting  (A2.19)  from  (A2.18).  For  X  ■  .6,  A  -  .4 
b<t+l)  -  b£c)  +  b(t)  +  .87  at+1 


b(t+l) 

m 


*&♦ 


[.01  (13  -  2m)  -  .4]  at+1,  m  4  12 

-.11  at+x ,  ra  ■  12 


Table  A2.1  gives  the  updating  coefficients  for  b  and  b  for  this  fitted 
&  r  °  om  m 


airline  model. 
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Table  A2.1 

Coefficient  of  a^+^  in  Updating  Formulae  for  Fitted  Airline  Model 


L 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

U 

* 

(t+i) 

\ 

om 

.58 

•  56 

.54 

.52 

.50 

.48 

.46 

.44 

.42 

.40 

.38 

.76 

.021/ 

(t+D 

> 

m 

-.29 

-.31 

-.33 

-.35 

-.37 

-.39 

-.41 

-.43 

-.45 

-.47 

-.49 

-.11 

.87^/ 

a/  Slope  b 


(t+1) 


b/  Intercept  b 


(t+1) 
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